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Evidence for orthorhombic distortions in the ordered state of ZnCr204: A magnetic 
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We present an elaborate electron-spin resonance study of the low-energy dynamics and magne- 
tization in the ordered phase of the magnetically frustrated spinel ZnCr204. We observe several 
resonance modes corresponding to different structural domains and found that the number of do- 
mains can be easily reduced by field-cooling the sample through the transition point. To describe the 
observed antiferromagnetic resonance spectra it is necessary to take into account an orthorhombic 
lattice distortion in addition to the earlier reported tetragonal distortion which both appear at the 
antiferromagnetic phase transition. 

PACS numbers: 75.25.-l-z, 75.50.Ec, 76.50.+g 



I. INTRODUCTION 

The intriguing physics of spinel compounds is in 
the focus of current solid-state research. The cur- 
rent hot debates on the origin of exotic phenomena 
and ground states in magnetic spinels concern, e.g., 
the Verwey transition in FeaO^ii^, heavy- fermion for- 
mation in 'LiY20^i^, colossal magnetoresistance in Cu 
doped FeCr2S4^i^, gigantic Kerr rotatioii^ and the or- 
bital glass state in FeCr2S4^, the spin-orbital liquid 
in FeSc2S4^, the colossal magneto-capacitive effect in 
CdCr2S4 and HgCr2S 4^°'^^ , the negative thermal expan- 
sion and strong spin-phonon coupling in ZnCr2Se4 and 
ZnCr2S 4^^'^'^i^^ , the spin dimerization in CuIr2S4i^ and 
yigT\20^, and the spin-Peierls-like transitions in 3- 
dimensional solid a^^'^^i^^ . The appearance of these fas- 
cinating ground states is attributed to the competition 
of charge, spin and orbital degrees of freedom, which are 
strongly coupled to the lattice. 

Additional complexity in the normal AB2X4 spinels 
arises from the frustration effects related to the topolog- 
ical constrains of the pyrochlore lattice of corner-sharing 
tetrahedra of the i3-site magnetic ions. In this geome- 
try the exchange interaction alone cannot select a unique 
ground state. As a result, the magnetic system remains 
in the disordered state down to temperatures much lower 
than the scale provided by the exchange interaction. In 
ZnCr2 04 strong direct antiferromagnetic Cr-Cr exchange 
is manifested by the Curie- Weiss temperature of about 
—400 K, while magnetic order appears only below 12 K 
via a first-order phase transition. At this transition the 
aforementioned degeneracy is lifted by a structural de- 
formation, which is reported to be tetragonalJ^ 

Neutron-scattering experiments have proven that non- 
collinear commensurate antiferromagnetic order is estab- 
lished below the transition temperature. However, the 
details of the magnetic structure are still under heavy 
debate. It was speculated that a multi-fc structure is 



formcdf^i^ Moreover, sample dependent intensities of 
the magnetic reflections suggest that ZnCr2 04 is criti- 
cally located close to several spin structures'^. 

In the present paper we report results of a magnetic 
resonance study in the ordered phase of ZnCr204. Ear- 
lier ESR studies'' concentrated on the paramagnetic 
state above T^ and did not report any resonance ab- 
sorption signals below T/v. We find several gapped res- 
onance modes, spin-reorientation transitions, and evi- 
dences for an orthorhombic structural deformation. We 
observe magnetic resonance signals originating from dif- 
ferent structural domains. We demonstrate that these 
domains can be effectively aligned by field cooling in a 
moderate magnetic field. The low-energy dynamics can 
be described qualitatively within the exchange-symmetry 
theory, indicating that magnetic order in ZnCr204 is gov- 
erned by a single order parameter. 



II. EXPERIMENTAL DETAILS 

ZnCr204 single crystals were grown by chemical trans- 
port reactions from polycrystalline starting material pre- 
pared by solid-state reactions of stoichiometric binary 
zinc and chromium oxides of 99.99% purity. Perfect sin- 
gle crystalline samples of octahedral shape and dimen- 
sions up to 3 mm on the edge were obtained. X-ray 
diffraction at room temperature revealed a single-phase 
material with the cubic spinel structure with a lattice 
constant a — 8.332(1) A and an oxygen fractional coor- 
dinate X — 0.263(1). The magnetic properties were stud- 
ied using a commercial SQUID magnetometer (Quantum 
Design MPMS-5) working at fields up to 50 kOe. 

Magnetic resonance measurements in the wide fre- 
quency range from 20 to 150 GHz were performed at 
the Kapitza Institute. For these measurements we have 
used a set of home-made transmission-type ESR spec- 
trometers equipped with a superconducting cryomagnet. 
High-sensitivity X-band (9.3 GHz) magnetic resonance 
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FIG. 1: (color online) Temperature dependences of the mag- 
netic susceptibility in different principal orientations. All 
curves are measured on cooling in the field of 10 kOe. 
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FIG. 2: (color online) Field dependences of the magnetization 
divided by field for the zero-field-cooled (ZFC, open symbols) 
sample and field-cooled (FC, closed symbols) for different ori- 
entations at T=2K. 



experiments were carried out using a Bruker "Elexsys 
E500 CW" spectrometer equipped with an Oxford In- 
struments helium gas- flow cryostat. Magnetic resonance 
absorption spectra were recorded at different frequen- 
cies for three principal orientations of the magnetic field: 
H||(001), (110), (111). The measurements were mostly 
done on zero-field cooled samples, the effect of field cool- 
ing was checked at certain frequencies. 



III. EXPERIMENTAL RESULTS. 

Figure [T] shows the temperature dependence of the 
magnetic susceptibility of ZnCr204 single crystals for the 
magnetic field applied along all three characteristic direc- 
tions of the cubic system. The susceptibility is isotropic 
above the antiferromagnetic (AFM) transition tempera- 
ture T/v = 12.5 K. At high temperatures the suscepti- 
bility follows a Curie- Weiss law, but deviates already at 
about 100 K and develops a broad maximum around 50 
K indicative for short-range AFM correlations. At Tn 
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FIG. 3: Field dependences of microwave absorption at differ- 
ent temperatures (ZFC sample). Inset: temperature depen- 
dence of the resonance fields. The narrow line at _ff = 13 kOe 
is a DPPH {g = 2.0) marker. 



one observes a discontinuous change of the data typical 
for a first-order transition. In the magnetically ordered 
regime the magnetic susceptibility of the zero-field-cooled 
(ZFC) sample shows a pronounced anisotropy with the 
highest value for measurements along the (001) direction 
(for the measurements in the field of 10 kOe). 

To investigate the anisotropy in more detail, the mag- 
netization was measured dependent on the magnetic field 
both for zero-field cooling (ZFC) as well as after field- 
cooling (FC) in 50 kOe. As shown in Fig. H the ZFC 
data manifest a non-linear behavior in small fields (up 
to 20 kOe) and a linear increase of the magnetization for 
higher fields. For H||(100) the M/H curve shows a linear 
increase indicating smooth rotation of the magnetization, 
while for the two other directions strong non-linearities 
are observed around 15 kOc, as typical for a spin-flop. 
Field cooling has only a weak effect for H||(001), but 
leads to nearly constant M/H for the other two orien- 
tations. In the field cooled sample the largest value of 
magnetic susceptibility is observed for H||(110). 

The evolution of the resonance-absorption spectrum 
with temperature is shown in Figure [31 At high temper- 
atures (in the paramagnetic phase) a single absorption 
component with a g-factor close to 2.0 is observed. The 
transition to the antiferromagnetically ordered state is 
clearly marked by the discontinuous transformation of 
the resonance absorption spectrum: below the Nccl tem- 
perature Tn the absorption spectrum consists of several 
components strongly shifted from the paramagnetic res- 
onance position. No hysteresis exceeding the resolution 
limit of 0.1 K was detected at the transition. On cool- 
ing below T^ the resonance lines first show a pronounced 
shift, but below 5 K the temperature dependence of the 
resonance positions is negligible. 

The shape of the resonance absorption spectra is 
strongly affected by field cooling. Figure [3] compares the 
resonance absorption measured on ZFC and FC samples. 
Here field cooling was performed at the field of 50 kOe 
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FIG. 4: (Color online) Comparison of the resonance absorp- 
tion observed at 1.8 K in zero-field cooled (ZFC) and field- 
cooled (FC) samples for different orientations of the magnetic 
field. Inset: Field cooling effect on the AFM soft mode at 
T = 4.2 K and / = 9.4 GHz. After FC only the background 
signal remains due to frozen oxygen). 



starting from 20 K. For H||(110) and (111) field cooling 
leads to the disappearance of some of the absorption com- 
ponents. The remaining absorption components are usu- 
ally slightly shifted from the corresponding absorption 
component measured on the ZFC sample. The vanish- 
ing absorption intensity does not necessarily add to the 
remaining components: for example, for H||(110) the in- 
tensity of the remaining component after field cooling is 
the same as for the ZFC sample. For H||(001) field cool- 
ing effects are less evident — all absorption components 
are observed in FC samples, field cooling leads only to a 
slight change of the absorption intensity. 

Actually, application of a field of 50 kOe during the 
cooling seems to be excessive. It is enough to cool the 
sample at the moderate field of 18 kOe to suppress some 
of the resonance modes as can be seen in the inset of 
Figure m Especially, the soft modes indicative for spin- 
reorientation, disappear after field cooling but reappear 
as the sample is rotated to another crystallographically 
equivalent position. Note that the nonlinearity of the 
magnetization curves vanishes after field cooling as well. 
The stability of the ZFC resonance absorption under pro- 
longed exposure to the magnetic field below Tn was also 
checked by the high-sensitive X-band measurements. At 
4 K the shape of the resonance absorption is reproducible 



to the finest details. However, at 8 K (which is still be- 
low Tjv) keeping the H||(110)-oriented sample at 18 kOe 
for 90 minutes leads to 30% reduction of the observed 
resonance absorption. 

The entire frequency-field diagrams for the different 
orientations of the magnetic field are given in Figure [5] 
These dependences demonstrate the presence of several 
resonance modes with zero-field gaps of 21 ± 2 GHz and 
113 ± 2 GHz. For H||(110) and (111) one of the res- 
onance modes softens in the magnetic field between 10 
and 15 kOc. Note that the nonlinearity of the magnetiza- 
tion curves is also observed at H = 10 ... 15 kOe for these 
directions (Figure [5]). Field cooling reduces the number 
of the observed resonance modes to two: one for each of 
the zero-field gaps. 



IV. DISCUSSION. 
A. Phase transition, domains, field cooling 

As documented above, the change of the resonance 
field at the phase transition is discontinuous. In conven- 
tional molecular field approximation, the shift of the an- 
tifcrromagnetic resonance field with respect to the para- 
magnetic resonance is proportional to the magnitude of 
the order parameter, i.e., to the sublattice magnetiza- 
tion. The discontinuous change of the resonance field at 
the phase transition indicates that the order parameter 
is not small even just below the transition temperature. 
This observation is in agreement with the first-order na- 
ture of the phase transition in ZnCr204 J^i^ 

The magnetic susceptibility of the paramagnetic phase 
is isotropic due to its cubic symmetry. Cubic symme- 
try is lost in the ordered state because of the lattice de- 
formation. The lattice-strain direction can take one of 
the equivalent crystallographic axes. The susceptibility 
tensor of the antiferromagnet is anisotropic, the orienta- 
tion of its principal axes is determined by the orientation 
of the order parameter fixed by anisotropic interactions 
with respect to the crystallographic axes. Therefore, the 
susceptibility tensors of different structural domains are 
oriented differently and the gain in the Zeeman energy 
is different for different domains. Thus, the application 
of magnetic field makes one of the domains more favor- 
able. It provides an obvious mechanism for the observed 
field-cooling effect and for the instability of the resonance 
absorption close to Tm^ described in the previous section. 
This assumption is in agreement with the increase of the 
magnetic susceptibility in the field-cooled sample (Figure 



HI). 

The formation of the domain structure in the case of a 
cubic-to-tetragonal lattice transition is well studied for 
ferroelastic systems. A complicated domain structure 
consisting of thin twinned domains is usually formed in 
ferroclastics (see, for example, recent Refs. [25li26l ). Twin- 
ning allows to avoid strong local strain at the contact of 
the domains with different directions of deformation axes. 



The thickness of twin domains observed in the doped fcr- 
roelastic HTSC compound YBa2Cu307 is about 10 — 100 
nm.— Twinning also leads to a slight tilting (of the or- 
der of Aa/a) of the domain axes from the correspond- 
ing crystallographic axes of the high-temperature phase. 
While the axes tilting is too small (Aa/a ~ 10^"^) to 
be of importance, small domain thickness could result in 
excitation of standing spin- waves with k ^ 1/L (i is a 
domain thickness) instead of a uniform fc = oscillation. 
This size effect can be a possible reason for the slight shift 
of the resonance absorption position in the FC sample. 
This allows to estimate the thickness of the crystallo- 
graphic domains. We assume a quadratic spectrum of 
antiferromagnetic spin- waves 



E 



A^ + a'^k'^ 



(1) 



here a ^ Ja. J denotes the exchange coupling con- 
stant between nearest-neighbour spins at distance a. If 
ak ^ Aq, the effective zero-field gap for the standing 
spin-waves is larger than for the uniform oscillation by 
S = a^k'^/{2Ao). Consequently, the resonance branches 
of the monodomain (FC) sample should be shifted down- 
wards (on the H-f plane) by S with respect to the res- 
onance branches of the multidomain (ZFC) sample: i.e. 
for the branches rising with the field the remaining ab- 
sorption component of the FC sample shifts at the given 
frequency to higher fields with respect to its position in 
the case of the ZFC sample. This slight shift is observed 
in the experiment (Figure S]), its magnitude does not ex- 
ceed the half-width of the absorption line. We estimate 
S/h as 1 GHz {h is the Planck constant). Then, the do- 
main thickness can be estimated as 
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(2) 



which yields (substituting J = 20 K,-2I Ao//i=20 GHz, ) 
i/a- 100> 1. 

In further discussion we assume that the domains are 
thick enough to be considered as bulk antiferromagnet, 
and that the domain walls do not contribute to the mag- 
netic resonance absorption. 



B. Application of the exchange-symmetry theory 

We use the theory of exchange symmetry^^ to analyze 
the experimental results. According to this theory, if the 
relativistic effects are much smaller than the exchange 
interaction, the order parameter can be represented by 
at most 3 unitary orthogonal vectors which transform 
by irreducible representations of the crystal symmetry 
group. The number of vectors and these representations 
define the exchange symmetry of the magnet. This ap- 
proach allows to describe all symmetry-based properties 
of a magnet without considering its detailed microscopic 
structure or any model assumption. For the case of a 



noncoUinear antiferromagnet the order parameter con- 
sists of at least two vectors. In case of only two vectors 
for the sake of simplicity we denote 1^'^-' = [l*-^-* x 1*^^^]. 

We derive the dynamic equations for homogeneous os- 
cillations from the Lagrange function 



^ = El(i^^^+T[l*'^xH] 



u„ 



(3) 



where 7 is the gyromagnetic ratio of the free electron, 
the constants li are related to the eigenvalues of the sus- 
ceptibility tensor, and the term Ua includes small rel- 
ativistic corrections to the main exchange part due to 
spin-orbital and dipole-dipolc effects. These corrections 
can be expanded by the components of the order param- 
eter. The first order of the Ua expansion is quadratic 
on If' . The second order terms can be omitted, because 
they are a^ times smaller, where a is the fine-structure 
constant. The Lagrange function must be invariant un- 
der all transformations of the crystallographic symmetry 
group, which results in some relations between the coeffi- 
cients in Ua expansion. These relations vary for different 
exchange-symmetry groups. 

The dynamic equations are obtained by taking a varia- 
tional derivative of the Lagrangian ^ over the rotations 
of the spin space. These equations should be linearized 
near the equilibrium position to obtain the eigenfrequen- 
cies of small oscillations. In general case, if the magnetic 
field direction deviates from the eigenvector of the sus- 
ceptibility tensor, the equilibrium positions cannot be de- 
termined analytically. To model the observed frequency- 
field dependences, we perform numerical calculations of 
the oscillation eigenfrcquencies. We use standard min- 
imization routines to find an equilibrium orientation of 
the order parameter. This modeling procedure is com- 
bined with a fitting algorithm using the constants li and 
the coefficients of the Ua expansion as fit parameters. 

When performing the expansion of the relativistic cor- 
rections, it is necessary to take into account that in the 
case of ZnCr204 the magnetic unit cell is larger than the 
crystallographic one»21i^. Therefore some components of 
the order parameter are not invariant under some of the 
translational elements of the crystallographic symmetry 
group. Since the crystal-symmetry group l?!^;, suggested 
in Ref . [2l|, has a point symmetry D2d in the vertex of the 
crystallographic cell, we will focus primarily on the point- 
symmetry subgroup. Note that this special property re- 
mains for all subgroups of -C'2d- Thus, on discussing the 
lowering of the lattice symmetry below T/v to D^, we will 
focus primarily on the point-symmetry subgroup. 

Although some representations of D2d allow weak fer- 
romagnetism, the susceptibility measurements do not re- 
veal any spontaneous magnetization. This can be either 
due to the spontaneous magnetization being too small, 
or, more likely, there is no weak ferromagnetism for the 
exchange group in the present case. Therefore, we will 
not take weak ferromagnetism into account in further 
discussion. 



C. Evidence for orthorhombic distortions below Tm 

Here we will demonstrate that the assumption of the 
tetragonal lattice symmetry in the ordered phase con- 
tradicts the experimental observation described in the 
previous section and the explanation of the experimen- 
tal findings requires a further reduction to orthorhombic 
symmetry. 

First, we note that the symmetry of the magnetic 
structure below T^ is lower than tetragonal. This state- 
ment follows directly from the observation of the distinct 
field-cooling effect for H||(lll), since this field orienta- 
tion is equivalent for all tetragonal domains. 

For the tetragonal lattice deformation only three types 
of different crystallographic domains, differing by the di- 
rection of the tetragonal axis z (z||(001), (010), (001)), 
can be formed at the transition. Since the magnetic sym- 
metry of the ordered phase is lower than the lattice sym- 
metry, two types of magnetic domains with different x 
and y axes can be formed in each crystallographic do- 
main. Here and further on in this paper we will denote 
by X and y the directions of the two-fold axes perpendic- 
ular to the z axis (5*4 axis for the D2d symmetry). Note 
that the reduction of the point symmetry from Oh to D2d 
could be done in two ways: (i) by removal of the [100] 
and [010] symmetry axes as well as (110) and (1 — 10) 
mirror planes; (ii) by removal of the [110] and [1 — 10] 
symmetry axes as well as (100) and (010) mirror planes 
with the former four-fold axes (100) and (010) becoming 
two-fold axes. In the second choice of axes the magnetic 
field aligned along the (111) direction would be equiv- 
alent for all domains and there would be no reason for 
the observed field-cooling effects. Thus, the first possibil- 
ity has to be realized. (IdentificationSi of the tetragonal 
phase as IAm2 also points to the first possibility.) There- 
fore, X and y axes are aligned along the diagonals of the 
cubic facets. 

Now all domains can be classified by the orientation of 
their "xyz" basis with respect to the cubic axes of the 
paramagnetic phase as shown in Table HI Some of these 
domains appear to be equivalent in the particular experi- 
mental conditions (here x, y, z are the unit vectors in the 
corresponding directions, trivial cases are combined): 



TABLE I: Classification of the AFM domains with respect to 
the cubic axes of the paramagnetic phase 



domain 


X 


y 


z 


(a) 


[110] 


[-110] 


[001] 


(b) 


[1-10] 


[110] 


[001] 


(c) 


[101] 


[10-1] 


[010] 


(d) 


[-101] 


[101] 


[010] 


(e) 


[Oil] 


[0-11] 


[100] 


(f) 


[01-1] 


[Oil] 


[100] 



H||[100] 

domains (e),(f): H||z 

domains (a), (b), (c), (d): H||x-t-y 
H||[110] 

domain (a): H||x 

domain (b): H||y 

domains (c),(f), (d), (e): H||^ + z 
H||[lll] 

domains (a),(c),(e): H||\/2x + z 

domains (b),(d),(f): H||V2y + z 

Note that for the H||(lll) field orientation there are 
only two types of different magnetic domains. As one 
can see from the experimental data in Figure [5] we ob- 
serve five resonance branches in this orientation: two 
originating from the higher gap, two originating from the 
lower gap and one in the high-field — low-frequency part 
of the frequency- field diagram. The observation of this 
fifth branch indicates the existence of a third zero-field 
gap which is not observed directly and is less or equal in 
magnitude than the lowest observed gap of 21 GHz. 

To apply the exchange-symmetry theory as described 
in the previous subsection it is necessary to specify the 
symmetry of the order parameter. The point symmetry 
group D2d has two-dimensional and one-dimensional ir- 
reducible representations. If l*-^-* and l*-^-* transform by 
a two-dimensional representation, then /i = I2 and the 
relativistic contribution to the Lagrangian ([3]) has the 
form 

+c(e)/f-/(i)/r)+iz?(/(3))' (4) 

If 1^^) and l*^^) transform differently under translations, 
then B = C = 0. In the other fundamental case, if l*-^-* 
and V'^' transform by one-dimensional representations, 
the relativistic contribution reads 



Ua = ^A(l^'^ 



2 1 



-B 



{e'y+c{iix'-ii''e')i^) 



or 



Ua = Ia (e')V Ib (eO' +^ (e^^f + li'^e^) (6) 

Depending on the coefficients of the Ua expansion, 
there are two general possibilities: (i) the 1*-*^ lie in the 
(xz) and {yz) planes; (ii) the l''^ lie in the planes in- 
cluding the z axis and the bisector of the (xy) plane (i.e. 
in the D2d group mirror planes). ZnCr204 corresponds 



to the first case, since in the second case for H||(lll) 
the order parameter vectors 1*^^^ l^^\ l*-^-* in all magnetic 
domains form the same angles with the field direction, 
which should result in the absence of the field cooling 
effects. 

We tried to fit the experimentally observed frequency- 
field dependences applying the exchange-symmetry the- 
ory as described above. However, we did not reach any 
reasonable agreement of the modeled curves with the ex- 
periment. The reason of this disagreement can be ex- 
plained qualitatively. Namely, it can be proven (see 
Appendix) that for the tetragonal D2d point symme- 
try the domain manifesting a spin reorientation transi- 
tion becomes indistinguishable from the other domains 
above the transition field. This is due to the fact that 
even for a magnetic symmetry lower than tetragonal, the 
anisotropic contribution Ua reflects the tetragonal crys- 
tal symmetry. Thus, after the spin reorientation tran- 
sition the reoriented domain appears to be in a state 
with the Zeeman energy lower than before reorientation, 
but with the same value of the anisotropic term Ua- The 
equivalence of the magnetic domains above the transition 
should provide the same magnetic resonance frequencies 
above the spin-flop transition. Moreover, there would be 
no reason for these domains to split again as the mag- 
netic fleld is reduced to zero. This, however, contradicts 
the experimental observation of the specific resonance 
modes corresponding to the domain undergoing a spin- 
reorientation above the spin-fiop field, as well as to the 
reproducibility of the low-field domain structure. 

The problem can be solved by the assumption that 
the lattice deformation at Tc involves not only a com- 
pression along the c-direction, but also a weak in-plane 
deformation leading to further reduction of the symme- 
try. There are two options: to exclude mirror planes 
(which yields D2 point symmetry) or to exclude second- 
order axes {C2v)- As we have shown above, the com- 
ponents of the order parameter lie in the {xz) and (jjz) 
planes. This limits our choice to D2 point symmetry 
(the corresponding space subgroup of IJjd is Dl (F222)) 
with second-order axes along (001) directions (z), and 
along (110) directions {x and y). This assumption pro- 
vides inequivalent crystallographic domains differing by 
the X and y directions. These orthorhombic distortions 
were not detected in the earlier structural studie o^^'^^ , 
most likely because they are probably smaller than the 
experimental resolution. 



D. Modeling of the AFM resonance frequency- field 
dependences for the case of orthorhombic distortions 

We modeled the antiferromagnetic resonance 
frequency-field dependences assuming orthorhombic 
distortions below Tjv. As was explained above, to write 
down the expansion of Ua we again will focus on the 
point-symmetry group. The D2 symmetry group has 
four one-dimensional representations. There are only 



four fundamentally different cases, the others can be 
reduced to them by renaming the axes x, y or vectors 

1(1), 1(2): 



Ua 



A 



C 



+E (/(D^f - /(D^)) + F (^(i)/(2) + /W^(2)) (7) 



Ua 



A 
2" 



B 

"2" 



-E(ii'¥J^^ 



/(i)/(2: 

i-z '■y 



A 



2 B 



C 



+ 



Z(i)/(2)) (8) 



(9) 



Here we again exclude weak ferromagnetism from the 
consideration. The first three cases are feasible only if l(i) 
and 1(2) transform in the same way under translations. 
The third case is feasible if l(i) and l(2) transform by the 
same one-dimensional irreducible representation. 

The low symmetry of the ordered state results in too 
many free parameters in the equations of spin dynam- 
ics (four to six coefficients in the Ua expansion and two 
of the three li constants). By fixing the zero- field gaps 
of the AFM resonance spectrum we can put only three 
constraints on these parameters. Other constraints are 
expected to appear during the fitting of the modeled 
AFM resonance spectra. This involves too many degrees 
of freedom for the assumptions on the sort of equilib- 
rium position, the way the spin-fiop transition takes, and 
the correspondence between structural domains and res- 
onance branches. An unequivocal quantitative reproduc- 
tion of the experimental curves was difficult to achieve, 
however, a reasonable agreement with the experimentally 
observed AFM resonance spectra was obtained for differ- 
ent sets of parameters and for the different forms of Ua- 

Here we present the results of modeling the case of Ua 
taken in the form ([8]). From Figure [5] one can see that the 
correspondence of model and experiment is fairly good. 
The correct quantity of resonance branches is obtained. 
The values of the zero-field gaps arc in good agreement 
with the experimental data. For orientations H||(110) 
and (111) the spin-reorientation transitions of certain do- 
mains arc well reproduced. 

Our modeling indicates that the frequency-field depen- 
dences can be reasonably reproduced assuming a single 
type of order parameter and taking into account all pos- 
sible domains. There are really just two types of dis- 
crepancies between the experimental data and modeled 
curves: first, there are predicted low-frequency modes at 
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FIG. 5: (color online) Frequency-field dependences of 
the resonance modes for the principal cubic orientations 
(open symbols-ZFC, closed symbols-FC). The numeric cal- 
culations (lines) are based on Eqn.(|8]). Fit parameters: 
7 = 2.8 GHz/kOe, 7i=3.931, 12=0.950, 13=1.00, A = 
65523.87 GHz^ B = 3444.56 GHz^ C = -6192.39 GHz^ 
D = -5664.49 GHz^ E = -8710.32 GHz^ F = 0. 



high fields which arc not observed; second, for the orien- 
tation H| 1(110) the modes corresponding to the domain 
with lowest energy (H||a;) are far from the experimen- 
tal data corresponding to the domain stable under field 
cooling (see the fit curves and experimental data marked 
by a, a', b, b' in the lower panel of Fig. [S]). These dis- 
crepancies can be due to the following reasons: at 9 GHz 
ESR experiments, fields up to only 18 kOe were available, 
whereas above 18 GHz the frequency-field dependences 
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FIG. 6: (color online). Upper frame: Calculated field de- 
pendences of the energy difference for various domains. The 
calculations are done for the same parameters as used for the 
AFM resonance spectrum shown in Figure [S] Middle frame: 
Calculated magnetization of the domains with unfavorable 
orientation. Inset: Calculated field dependences of the Euler 
angles for the H||(110) y-domain. Lower frame: Calculated 
magnetization of the domains with favorable orientation of 
the order parameter. 



of the low-frequency modes in question are rather flat. In 
the case of I2 = I3 they become field independent which 
results in a strong broadening of the resonance modes in 
the experiments at fixed frequency. Thus, these modes 
cannot be detected in a field scan. Concerning the dis- 
crepancy for the H||x domain, note that the modeled 
curves still demonstrate an " anti-crossing" which can be 
enhanced, for example, by higher order terms in the en- 
ergy expansion. 

Our modelling demonstrates that the spin- 
reorientation transition can be realized by a continuous 



rotation of the order parameter as shown e.g. for the 
y domain in the inset of Fig. [HI This explains the 
absence of a sharp change of the magnetization at 
the spin-reorientation transition which is simulated in 
the middle and lower frames of this figure. The field 
dependence of the magnetization is qualitatively well 
reproduced: it demonstrates a nonlinear behavior with 
a characteristic change during the spin-reorientation 
for the unfavorable domains, which would disappear 
after field cooling, and an almost linear behavior for 
the favorable domains which remains after field cooling. 
The energy difference between different domains shown 
in the upper frame of this figure explains the observed 
field-cooling phenomena. The domains corresponding to 
the modes surviving the field cooling have the lowest en- 
ergy in the magnetic field. For the orientation H||(lll) 
the calculated energy difference between the domains is 
similar and, thus, the field-cooling effect is comparable 
to the case H||(110). For H||(001) this difference is much 
smaller which can probably explain the much weaker 
field-cooling effect for this direction. However we have 
to admit the problem in reproducing the experimentally 
observed anisotropy of the field-cooled magnetization 
at low fields: the experimentally highest susceptibility 
for the monodomain FC sample is observed for H||(110) 
(Figure [2]), while the modelled curves for the presented 
parameter set indicate a slightly higher susceptibility for 
the most favorable domain in the H||(100) (Figure [6]). 
Note that the modeled difference of the three favorable 
domains is tiny, thus for low fields some weak effects 
related to shape anisotropy or surface effect may distroy 
the stability of the uniform magnetization. 

Summarizing the results of our modeling, we conclude 
that the low-energy dynamics of the ordered phase of 
ZnCr204 can be described by the assumption of a single 
type of magnetic order below T^y. This puts a question 
mark on the possibility to realize different spin structures 
in ZnCr204 as suggested earlier—. We suppose that the 
complications in the determination of the magnetic struc- 
ture of ZnCr204 by neutron scattering were caused by 
the unaccounted orthorhombic deformations and effects 
of multiple domains. 



V. CONCLUSIONS 

We have performed a detailed study of the low-energy 
dynamics of the ordered phase of the frustrated antifer- 
romagnetic spinel ZnCr204. We have proven directly 
that multiple domains exist below the transition tem- 
perature Tn- We have demonstrated that some of these 
domains are effectively suppressed by field cooling. Spin- 
reorientation transitions indicated by softening of certain 
antiferromagnetic resonance modes and by the nonlinear 
behavior of the magnetization are observed. 

These results are incompatible with the earlier pro- 
posed symmetry of the distorted lattice. Thus, we con- 
clude that the lattice deformation at the phase transition 



involves small in-plane distortions besides the tetragonal 
distortions. We suggest that the actual symmetry of the 
lattice below Tm corresponds to the orthorhombic D^ 
symmetry. 

We have demonstrated that the low-energy dynam- 
ics can be reasonably described within the framework of 
the exchange-symmetry theory, assuming a noncoUinear 
magnetic ordering characterized by a single order param- 
eter. 
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APPENDIX A: INSTABILITY OF THE CERTAIN 

DOMAINS IN CASE OF THE TETRAGONAL 

LATTICE SYMMETRY 

Here we will prove that for the D2d point symmetry 
at the spin-reorientation transition the domain demon- 
strating this transition becomes indistinguishable from 
the other domains. 

First we consider the case when 1^^^ and 1^^^ transform 
by two-dimensional representation of D2d- Minimizing 
the energy @ we get the orientation of the order param- 
eter in the absence of magnetic field. For \B\ > \A\ the 
solution consists vectors lying in the mirror planes, which 
is not the case of ZnCr204 (see above). For \A\ > \B\ 
and A < the two solutions are I'-^^Hx, l^^^ lying in the 
(yz) plane and l^^-'Hy, l'^' lying in the {xz) plane. The 
case oi A > can be reduced to this one by renaming the 
axes X and y. These solutions define two different mag- 
netic domains. If the susceptibility along l^'^^ is less than 
along l'^) and l'^', then for -ff||x one of the domains is 
already in its minimum of the Zeeman energy. The other 
one is not in the minimum and at some value of magnetic 
field it will undergo a spin-reorientation transition. As 
it turns out, there is only one possible state for it after 
the transition, the same as for the first domain. So after 
the spin-flop these two domains will be indistinguishable. 
If the susceptibility along l'^^ and l'^^ is less than along 
l'^', both domains are not in the minimum of the Zeeman 
energy when the field is applied along the x axis. Still 
after the spin-flop both domains again become indistin- 
guishable. 

Now we consider case when l'-'^' and 1^^' transform by 



onc-dimcnsional representations of D2d- To define the 
orientation of the order parameter we minimize the en- 
ergy ([5]). We get the result that one of the vectors (let 
it be l'^^) is aligned in the {xy) plane, but its orienta- 
tion in this plane remains arbitrary. To find a solution, 
it is necessary to take into account the next order terms 
in the Ua expansion. There is no need to write down 
all of them, just note that due to the tetragonal sym- 
metry the dependence of Ua on the angle (f> between 1^^^ 
and X is F cos (40). There are two sets of solutions de- 
pending on the sign of F. For F > the solutions are 
(j) = 7r/4, 37r/4, 57r/4, 77r/4, i.e. l'^^^ lies in the mirror 
plane, which is not the case of ZnCr204. For F < 
the solutions are </> = 0, 7r/2, tt, 3n/2. These solutions de- 



fine two magnetic domains: for one of them l^^^ | |x and 
for the other l*-^-'||y. First, no matter along which l''' 
the susceptibility is largest, a spin-reorientation transi- 
tion is expected for H||(x ± y) (i.e. H||(001)). This 
spin-reorientation is either rotation of the order param- 
eter around the z axis by 7r/4 or, in the special case of 
the largest susceptibility being along the z axis, rotation 
of the largest susceptibility direction to the (xy) plane. 
However, such a transition is not observed in our ex- 
periments. Second, the spin-flop transition observed at 
H| 1(111) can be caused only by rotation of the order pa- 
rameter around the z axis, but after this rotation both 
domains become indistinguishable. 
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